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Abstract 

We consider the production of high transverse momentum gluons in the McLerran- 
Venugopalan model of nuclear structure functions. We explicitly compute the high 
momentum component in this model. We compute the nuclear target size A depen- 
dence of the distribution of produced gluons. 
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1 Introduction 



Understanding parton distributions formed at the initial stages of a collision between 
two heavy ions is a very important open problem. In a previous paper we have set 
up a formalism for the calculation of the gluon distribution in the framework of 
the McLerran-Venugopalan model of the nuclear structure functions |l],|2]. In this 
approach the valence quarks in the nuclei are considered as classical sources of color 
charge. This picture is similar to that developed by Mueller for the gluon structure 
functions for heavy quark systems 0. The "initial values" of the color distribution 
of the valence quarks in the two colliding nuclei is given by 

Jl 2 (z(x)) = 6 v± d(x*)g Pl , 2 (x ± ) (1) 

As is apparent from the the light cone delta functions S(x Zf ), the quarks in the nuclei 
appear as infinitely thin sheets of nuclear matter moving at the speed of light in 
positive and negative z directions respectively. 

These color charge distributions generate a classical glue field according to 
classical Yang-Mills equations 

[D„,F^] = Yl U[A](x,z m (x))J^(z m (x))U[A](z m (x),x) (2) 
m=l,2 

Here zi^ix) = x\ x ± =0 serves as a reference point used to define the initial value of 
the charge distribution. Due to the covariant current conservation [D u , J u (x)] = 
this initial distribution evolves along the trajectory of a particle via parallel trans- 
port. This is the origin of the link operators 

U[A](x,z m (x)) :=Pexp-ig T dw^{u m ) (3) 

J Zm(x) 

connecting the initial point z m (x) and the point x on the trajectory uj m , which 
appear on the right hand side of eq. (Q) . 
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The gluon distribution function is denned in terms of classical solutions of eq. 
@, and is related to the following quantity (for precise definition see Section ||) 

< Mk)A*(k) > p (4) 

Here A(k) is the Fourier transform of the classical solution. The averaging over the 
color charge distributions is performed independently for each nucleus with equal 
gaussian weights 

(0) p = J d Pl dp 2 O exp j d 2 x ± tr[pl(x L ) + f%(x±)]} (5) 

The purpose of this paper is to calculate this distribution function perturba- 
tively, to lowest nontrivial order in the inverse powers of the transverse momentum 

Ctsfl/k. 

Basic input in our approach are the classical fields generated by single nuclei 
which have been derived earlier in [|], f|,|5|| . These solutions remain valid before the 
collision and yield initial conditions for the field after the collision. The solution in 
the one nucleus case is of the formf] 

Am = 

Al = 9(x*)a\ >2 (x ± ) (6) 

The functions oti are implicitly determined by the "dimensionally reduced" ver- 
sion of the Yang- Mills equations 

a4 = -^U m (x ± )d l Ul(x ± ) 
d l a l m = gp m (x±) (7) 

An obvious property of the solutions eq.@ is that the "transverse" components of 

the field strength vanish F 13 = 0. The gauge potentials themselves vanish in front 
1 On these classical solutions the link operators in eq. (^) drop out. They may however be 
important if one starts to consider quantum corrections in powers of a s . 
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of the moving charge and are a pure gauge behind it. The only physical information 
is contained in the discontinuity at the worldline of the quarks which generate the 
gluon fields. The field strength does not vanish only on these worldlines and is 
confined to infinitely thin sheets. 

Let us now turn to nucleus-nucleus collisions. Obviously the single nucleus 
solutions are still valid everywhere except in those regions of space-time which are 
in causal contact with the collision point, i.e. in its forward light cone (see Fig. |l]). 
In the abelian case even there the solutions would be directly given in terms of the 




Fig. 1: Regions with different 
structures of the gauge poten- 
tial: 

In regions 1 and 2 we have the 
well known one nucleus solu- 
tions 0:1,2- While the gauge 
potential in the backward light 
cone is vanishing we have a non- 
trivial solution in the forward 
lightcone, region 3 



single nucleus solutions as a sum of two pure gauge fields which again constitutes 
a pure gauge field. Hence in this case there would be no nontrivial effects on the 
classical field level. To generate photons and field strength one would have to 
resort to quantum effects. This is fundamentally different in the non-Abelian case. 
Obviously now, the sum of two pure gauge field is no longer a pure gauge field 
due to nonlinear effects, and therefore it does not solve the Yang Mills equations 
for the two nuclei problem. Nevertheless the one nucleus solutions provide initial 
conditions for the evolution of the gauge field in the forward light cone as discussed 



4 




In the remainder of this paper we will adopt the gauge condition 

x + A~ +x~A + = (8) 
In this gauge the fields in the forward light cone can be cast in the form 

A^ 1 = ±x a(r, x±) 

A 1 = ai(r,x ± ) (9) 

where r = V2x + x~ is the proper time. Note that this representation does not 
introduce any physical assumption not already present in the single nucleus solution. 
The relation between is a consequence of the gauge condition. Further, the fact 
that the functions a(r, x±) and (Xj_(t, x±) do not depend on the space time rapidity 
T) = l/21nx + jx~ is a natural consequence of the absence of a longitudinal length 
scale in our initial conditions with its infinitely thin nuclei moving at the speed of 
light. 

Since the valence color charge densities vanish inside the light cone, the fields 
A^ of eq. @ should solve the homogeneous equations 

[D^F^] = (10) 

The initial conditions are determined by patching the solutions in the three regions 
(see Fig. [l]) together and eliminating the discontinuities at r = 0. This specifies 
the forward light cone fields a^_| T=0 and a| T=0 in terms of the known single nucleus 
solutions outside the forward light cone. 
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2 Perturbation theory 



Using the form (g) of the gauge potential, the Yang-Mills equations ( |10| ) in the 
forward light cone are 

1 



- t d T T d T a + 



A, 







1 



Di,d T a\_ +igT[a,d T a] = 



—d T Td T a\ — igr 



(11) 



Matching conditions on the solution lead to the initial conditions for a(r, and 
aj_(r, in terms of the single nucleus solutions a l m : 



a 



r=0 



T=0 



a| + a 2 



!£ 
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(12) 



In this paper we construct solutions in the weak field limit by expanding first 
the initial conditions and then the fields within the forward light cone in powers of 
P- 

Let us first concentrate on the initial conditions. To determine the initial condi- 
tions perturbatively, we need the single nucleus fields ot % m , i = 1, 2 to the appropriate 
order. 

Recall that for a single nucleus, the Yang-Mills equations reduce to the two 
equations (H). The first states that a % m is a pure gauge field. The second equation 
determines a l m in terms of the charge density p m (x±) of the individual nuclei via 



a. 



ig 



gpm{x±) 



(13) 



To second order in p, the solution of ( |i3[ ) is 



0( P 3 j} (14) 



where we have defined 

± 9 



Pm (15) 



Vl 2 ' 

Note that the first order term in (|l4] ) is longitudinal, whereas all higher corrections 
are transverse. Their sole purpose is to render a l m a pure gauge field. Plugging 
( |i~H ) into (12) gives us initial conditions up to second order in p m . 



To solve the equations of motion ([ll]) in the forward light cone perturbatively, 
we also have to expand the gauge fields there in powers of p m 



a = a (n) 

n=0 

oo 

«i = E«i(n) (is) 



n=0 

First order: 

To first order the equations of motion are necessarily linear 

^3<9 T T 3 <9 r a (1 ) - V_L 2 a ( i) = 
d^ra^ = 
UrTdra^ - (Vi 2 ^' - &&) aj {1) = (17) 

They are to be solved subject to the initial conditions 

<X(i)(t,x±) | t=0 = 
«i(i)(r,a;±) | r=0 = - & (<? + <f) (x ± ) (18) 

The initial condition for the "±" -components to this order is trivial. 

With these initial conditions, the first order solutions are completely specified 
and essentially trivial. We obtain vanishing am 

«(i)(r,x ± ) = (19) 



and r-independent ctj_(i) 

ai ( i)(r,x ± ) = -^(V + <£ 2 )(* ± ) (20) 

which obviously is pure gauge. 

To this order the problem is structurally identical to its abelian counterpart. 
Commutators do not appear neither in the equations of motion nor in the initial 
conditions (these contributions are necessarily at least second order in p m ), and 
hence no field strength is generated in the forward light cone. 



Second order: 

In the second order this changes immediately. Using the first order information 
almost all nonlinear terms drop from the second order equations which now read 

d T a {2 ) - Vi_ 2 a( 2 ) = 

&d T ai_\ 2) = 

^d T rd T a ± \ 2) - (Va 2 ^ - a ± l 2) + ig ^[aj^a^] = (21) 

The only inhomogeneity in these equations comes from the - pure gauge - a _L(i)- 
In fact using the residual (r-independent) gauge freedom, we can remove the inho- 
mogeneous piece altogether. 

We define by 



A„ = V(x ± ) 



ig 



V-\ X1 _) (22) 



where V(x±) = exp — igcf)±(x±). </>j_(x_i_) is to be determined such that it removes 
nonlinear terms from the equations of motion for e at least to second order. In fact, 
it turns out, that it is possible to remove the inhomogeneous terms up to third order 
in p requiring 

0*^=0 = (23) 



In higher orders it is no longer possible to linearize the equations. Using equations 



(Eg) and (23) we find for 



+ (f>2 + 0(p 2 ) 



With this choice of V the first order solution vanishes 



(24) 



e (i) 



= 
= 



(25) 



The second order equations become 

1 



d r T 3 <9 r e (2) - V ± 2 e ( 2) 



d l d T e 



-d T Td T e\ 2) - (s ij V± 2 - e 



(2) 

3 

(2) 



The initial conditions on a 1 , and a are transformed into 









-(2) k=0 



-i<l [SV-d^d* ) <. n .<>>o 2 



(26) 



(27) 



\2) |t=0 



!£ 

2 



(28) 



Obviously, the gauge condition 3* e* | T= o = in this order is preserved for all proper 
times r. Hence the "i" -components of the gauge field may be represented as 



(29) 



Using e and Xi t ne second order Yang- Mills equations and initial conditions simplify 
even further. In terms of these variables, we finally obtain 



1 



1 



d T r 3 d T e {2) - V ± 2 e 



£ (2) 



d T TdrX(2) ~ Vl 2 X(2) 








(30) 



and 



£(2) \t=0 
X(2) |r=0 



2 



&<h,&<h 



-%g e 



(31) 



We stress that, even though the equations of motion ( pOD linearize to this order 
due to our gauge choice, the intrinsic nonlinearities of the system are very impor- 
tant. They show up in the initial conditions (|3l] ) which carry nontrivial physical 
information. 



3 Perturbative solutions 



To solve (|30|) , we note that the equations are of form 
4-^r m r /(r,^)-V ± 2 /(r,^) 



(32) 



where in our case m = 1,3, / = x-, e - Factorizing r and x^-dependence, we imme- 
diately identify the eigenfunctions 



(33) 



where ui^ = Jk\ and Z u may in general be any linear combination of the Bessel 
functions J U ,N V . 



The initial conditions (31) force our solutions to be regular at r = and thus 
preclude any admixture of Neumann- functions N v . 



A general solution is therefore of the form 
f d 2 k± 

F m (r,x±) = jj—^h m (k±)f kx (T,x±) 

s lfc i^-^ J h m (y±) M 1 -^ 1 -W,(u;t) (34) 



■ d 2 k A 
d 2 k±d 2 y_i 
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For equations (|30|) this yields 



X(2,(t,*l) = e«-»1 k,( V1 )AW (35) 

with /i3 and /ii determined by the initial conditions eq.([H|) as 



hi(x±) = -ig=~5 d l 4>i,9>4)2 (x±) (36) 
Vl 



In the following we will need the large r asymptotics of the solutions: 

X(3)(r,x ± ) = /^e^M*i) {^cos(o,r-?[) + 0(^)1(37) 



e {2) (r,x ± ) = I^L e ^h 3 (k ± ) cos(o,r - ^) + 0(^)1 (38) 



4 Distribution functions 

The solutions of the Yang Mills equations at asymptotically large r are of the form 

e °'* (T ' a: - L) = J0? K '^{ ailk ^ e ' tl ^~'" T + aa } (39) 

Here u> =\ k±\, and k 1 = e^k 3 /to. The notation C.C. means complex conjugate. 

To derive an expression for the energy density, we recall that r is large. Near 
z = 0, this implies that in the range of z where r ~ t 3> z the solutions are z 
independent. This means they asymptotically have zero p z . Now suppose we are at 
any value of z, and r is large but t ~ z. We can do a longitudinal boost to z = 
without changing the solution. Again in this frame the solution has zero p z . We 
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see therefore that for the asymptotic solutions the space time rapidity is one to one 
correlated with the momentum space rapidity, that is at asymptotic times we find 
that 

V = ^M x+ / X ~) = y = ^ln(p + /p~) (40) 

To proceed further, we compute the energy density in the neighborhood of 
z = 0. Here asymptotically r = t. The energy in a box of size R in the transverse 
direction and dz in the longitudinal direction, with L < t becomes || 

dE = d i /£^£i«fe)i 2 ( 41 ) 

^ ' i,b 

Recalling that dy = dz/t, we find that 

dE 1 



dyd 2 k ± (2vr) 3 



uJ2\^(k ± )\ 2 (42) 



and the multiplicity distribution of gluons is 

dN 1 dE 



(43) 



dyd 2 kj_ ui dyd 2 k± 

As we expect for a boost covariant solution, the multiplicity distribution is rapidity 
invariant. 

This is now easily compared to the asymptotic behavior of the second order 
solutions pf), We read off that 



1 h 3 (k ± ) 



IT UJ 



aa(fci) = -7=«w/ii(fc±) ( 44 ) 

V7T 



Hence we get the distribution function 
dN 1 



dyd 2 k± uj dyd 2 k 

1 



(27T) 



2 tr ai(kj 



3 







2 f/i 3 (fc±)4(fc±) 



£ _ tr ^ -U + ^(fe^T^) (45) 



(27r) 3 . 7T UJ 
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Now we use that 



h 3 (kj 
hi(k_i 



d y±e "'"- L a± - ig ■ 



{[avi,av 2 ]}( y± ) 



d y±e~ 



J-L 



(46) 



and get 



dyd 2 k\ 



9 2 2 
(2vr) 3 vr*; 2 



'2„, „ifc i -x i 



tr 



aVi,^> 2 l (s±) [9*01,^1 (0) 



(47) 



where, Sj_ is the transverse area of the system. As indicated, we have to average 
over p with the Gaussian weig ht eq.( (§)). Using (gj) 

we find 



, , , I'fAji P i P±-( x ~v)± 



(2tt) 2 V \ 



(48) 



and the above turns into 



dN 



dyd 2 ki 



fiijgkl , e ij e kl 
d^ 

(2vr) 2 p^p + kf 



X (2^)3^ 

4, /r a _6ir_& _ai\ [ d*p ± pY (p + W ( P + k) 1 



<? 6 2 



•/iV c (iV c 2 - 1) 



1 



(49) 



x (2tt) 3 vr*; 2 r ' ' c v ' c ^7 (2tt) 2 p 2 (p + k) 2 

The integral in equation (^) is infrared divergent. This however is an artifact of 
the weak field expansion employed above. As discussed in Jl|,§,§, ||, the inclusion 
of higher order terms for the initial fields would generate a mass scale of order a s fi 
and hence regulate this divergence. With logarithmic accuracy we obtain therefore 



dN 



dyd 2 k L x (2vr) 3 vr 

This is the main result of this work. 



S ±7 ^ViV c (iV c 2 -l) 



1 1 



hi 



A- 



2tt ki {a s pY 



(50) 
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5 Conclusion 



In this paper we have calculated the high k± asymptotics of the distribution func- 
tion of gluons produced in a collision of two ultrarelativistic heavy nuclei in the 
framework of the McLerran-Venugopalan model. 

Our result, eq. (|50|) has a very natural interpretation in analogy to the naive 
parton model. The number of produced gluons in the parton model is given by 

N(k±) = J ^p±N 1 (p ± )N 2 (k± - p±)a( P± , k ± ) (51) 

where N m (q±) is the probability to find a gluon with transverse momentum qj_ 
in the m-th nucleus, and a is the cross section for production of the gluon with 
momentum k± in the final state in the collision of the two gluons. To the lowest 
order in perturbation theory the cross section is 

cv 

a(k±) oc ^ (52) 



The distribution functions for the one nucleus case N m were calculated in |1|,||] and 
were found to have the same functional form as in the lowest order in perturbation 
theory 

N m (q ± ) oc ^ (53) 
1± 



Substituting this into eq.(|5l|) gives the result of eq.(|50|). 

Indeed the calculation performed in this paper can be represented in terms of 
Feynman diagrams. The gluon distribution function of eq.(|50|) has a representation 
of Fig. | 

The reader may have noticed one peculiarity in the preceding discussion. Al- 
though we have been using the parton model notations, some of the gluons involved 
in the process depicted on Fig. @ are virtual. To produce the real final state gluon, 
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Fig. 2: Diagrammatic rep- 
resentation of the lowest order 
distribution function 



the internal gluons must be off shell with virtuality of order of their momentum. 
In fact if one wants to think about this process in terms of real partons, it must 
be represented as "two quarks go into two quarks and a gluon", rather than "two 
gluons go into one gluon" , which is kinematically forbidden for on shell particles. 

As we have seen, to this order our calculation essentially reproduces perturba- 
tion theory. The reason is that we have expanded in powers of the valence charge 
density, and thereby have ignored the fact that the classical fields involved are ex- 
pected to be strong. The real nonperturbative nature of this approach will become 
apparent when those are fully taken into account. Some of the qualitative con- 
sequences of such an analysis however, can be understood already at this point. 
Consider for example, the total multiplicity of produced gluons. 

"-I^-Jk < 54 > 

Since dy d j2 k± is proportional to the transverse area of the system, and is a function 
of the dimensionless ratio /J?/k\, the integral has to be proportional to Sj_/U 2 . Re- 
membering that for large nuclei \i scales with the number of nucleon in the nucleus 
A as fj? oc A 1 / 3 , we conclude that for large nuclei the total multiplicity scales as 

TV oc A (55) 

On the other hand, the total multiplicity at momenta larger than some fixed mo- 
mentum p 3> an will coincides with the perturbative result 

N = J d 2 k L ^-ocS ±t i 4 ocA 4 ^ (56) 
15 




Therefore it is very important to go beyond the weak field expansion and in- 
corporate truly nonperturbative effects due to strong classical fields. 

In comparing our results with the parton cascade model of Geiger et al 0-^] one 
should keep in mind the following differences between the two approaches. In the 
parton cascade, the leading order in a s scattering occurs by gluon gluon scattering. 
This populates the high p± tail of the distribution by scatterings of low p± gluons. 
In our case, we have an intrinsic p± for the gluons and the gluons are far enough 
off mass shell so that we can produce high pj_ gluons by glue-glue goes to single 
gluon scattering. Consequently, although the p± dependence of the two results is 
similar, in the high p± region we have a different dependence on as- We have one 
power less, since the factor a 2 /j? in eq. ( |50| ) just sets the scale of the intrinsic single 
nucleus glue distribution. 

The scattering processes, which are leading in our approximation are also 
present in cascade codes as composite processes qq — ► qqg, where some of the 
quarks are off shell by an amount p\. Although these contributions are formally 
of order a 3 , for very large nuclei the quark distribution functions will be large and 
these processes may become leading. This seems to be the natural point of contact 
between the two approaches. 

Of course at some p± 3> /x, our approximations for computing the gluon dis- 
tribution function break down. This presumably first occurs when gluon brems- 
strahlung softens the high p± distribution from its l/p± behavior, and the spec- 
trum steepens. The main point is however not what happens in the tail of the 
distribution. This is simply where it is most easy to compute. In the center of the 
p± distribution of produced gluons, the gluons arise from a non-linear evolution of 
the gluons fields. This is cause by the quantum mechanical nature of the initial 
state and the charge coherence of the initial state interactions. In this case, the 
final state distribution of gluons bears scant resemblance to that of the initial dis- 
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tribution convoluted over hard gluon scattering. Before one knows whether there 
are substantial quantitative differences, one must of course numerically solve the 
problem for gluons in the center of the distribution. 
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